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can be infrared divergent — but that they are infinitely suppressed relative 
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1. Introduction 

The path integral in general relativity is a sum over geometries, and it is natural to 
ask whether this sum should be extended to include topologies as well 0. This question 
can take two forms: (1) should path integrals include sums over intermediate spacetime 
topologies ("spacetime foam"); and (2) should we allow transitions between different spatial 
topologies? These issues are closely related — it is hard to imagine a formalism that permits 
spatial topology change while forbidding sums over intermediate topologies — but they are 
distinct. In this paper, we focus on the second. 

Topology change, if it occurs, is expected to be a quantum mechanical process 0. Un- 
fortunately, we do not yet have a viable quantum theory of gravity with which to compute 
topology- changing amplitudes. We can, however, look for models that give us hints about 
what to expect in the full theory. 

One particularly interesting model is (2+l)-dimensional gravity, standard general rela- 
tivity in two spatial dimensions plus time. The underlying conceptual issues of quantum 
gravity, and some of the technical aspects as well, are identical in 2+1 and 3+1 dimensions. 
But the elimination of one dimension vastly simplifies the theory, making many computations 
possible for the first time. Moreover, general relativity in 2+1 dimensions is renormalizable 
(it is in fact finite), allowing us to avoid the difficult problems of interpreting path integrals 
in (3+l)-dimensional gravity. 

Indeed, as first noted by Witten ||,f§], the path integral for (2+l)-dimensional gravity can 
be evaluated explicitly in terms of a standard topological invariant, the Ray-Singer torsion. 
In first-order form, the action may be written schematically as 



(see the next section for details). The Lagrangian is cubic in the fields, but transition am- 
plitudes can still be computed nonperturbatively: the triad e acts as a Lagrange multiplier, 
giving a delta functional that absorbs the remaining integral over uo. With careful gauge- 
fixing, one obtains a combination of determinants of Laplacians known as the Ray-Singer 
torsion, a well-known topological invariant of a flat bundle over M. This torsion, in turn, 
is equivalent to the Reidemeister torsion, a combinatorial invariant that can be relatively 
simple to compute. 

In reference [|J, Witten points out that this equality makes the explicit computation of 
topology- changing amplitudes possible. He obtains a set of selection rules relating holonomies 
of u on the initial and final boundaries of M. In reference [[|, Amano and Higuchi obtain 
a much more stringent set of selection rules by demanding that the initial and final bound- 
aries be spacelike. They further write down a set of three-manifolds, expressed in terms of 
Heegaard decompositions, that interpolate between all topology changes allowed by their 
selection rules among surfaces of genus greater than one. 

The initial goal of this paper was to compute the Reidemeister torsion for the simplest of 
these interpolating three-manifolds, thus obtaining an explicit topology-changing amplitude. 
We do so in section 6. In the course of the computation, however, we found it necessary to 
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reanalyze Witten's expression for the path integral, treating boundary conditions somewhat 
more carefully. To our surprise, we discovered that the path integrals for the Amano-Higuchi 
manifolds typically involve zero-modes of the triad e, leading to infrared divergences of the 
type discussed by Witten in "classical" spacetimes [§]. 

To make sense of our topology-changing amplitudes, we must therefore add an infrared 
cutoff and compare the results to similar topology-preserving amplitudes. Unfortunately, the 
meaning of the word "similar" is not entirely clear. As noted above, once we permit spatial 
topology change, we ought to also allow sums over intermediate spacetime topologies as well. 
Such sums are typically badly divergent |J, and it is not clear how to normalize the total 
amplitude. This is a deep conceptual question, which goes beyond the issue of divergences: 
a topology-changing amplitude relates states in different Hilbert spaces, and we do not yet 
have a "Fock space" that includes multiple topologies. 

We shall not attempt to address this issue here. Instead, we take the less ambitious 
route of normalizing our amplitude relative to a topology-preserving amplitude coming from 
a three-manifold with a related Heegaard decomposition. We shall show that the effect 
of the infrared divergences discussed above is to infinitely suppress the topology-changing 
amplitude. In the absence of a clear prescription for normalization, we cannot yet claim that 
quantum gravity prohibits topology change, but we find this result rather suggestive. 



2. Action, States, and Boundary Conditions 

We begin with a brief review of (2+l)-dimensional gravity in the first-order formalism B. 
Let M be a three-dimensional spacetime. As our fundamental variables, we choose a triad 
— a section of the bundle of orthonormal frames on M — and an SO (2,1) connection on 
the same bundle, which we describe by a connection one- form uo^ a b . The standard Einstein- 
Hilbert action is 

/ g rav[M] = e a A (duo a + -e abc uj b A u c ^j , (2.1) 

where e a = e^ a dx^ and u a = ^e^u^cdx^. The equations of motion coming from this action 
are easily derived: they are 

de a + e abc u b A e c = (2.2) 

and 

R a = du a + ^e abc u b Aw c = 0. (2.3) 

The first of these is the condition of vanishing torsion, which ensures that the connection u 
is compatible with the metric g^ v = ri ab e il a ej' . The second is then equivalent to the ordinary 
vacuum field equations of general relativity. Note that ( |2.3| ) can also be interpreted as the 
requirement that the connection uj be flat, a feature peculiar to 2+1 dimensions that goes a 
long way towards explaining the model's simplicity. 

The action (|2.1|) is invariant under local SO (2,1) transformations, 

Se a = e abc e b r c 

5u a = dr a + e abc uj b T c , (2.4) 
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as well as "local translations," 

Se a = dp a + e abc u b p c 

5uj a = 0. (2.5) 

These transformations together form an ISO(2,l) algebra, and Witten has shown that the 
one-forms e a and uo a constitute an ISO (2,1) connection. J grav is also invariant under dif- 
feomorphisms, but this is not an independent symmetry: when the triad e^ a is invertible, 
diffeomorphisms in the connected component of the identity are equivalent to transforma- 
tions of the form (|2.4|) - (|2~5|) . We must still account for the "large" diffeomorphisms — the 
mapping class group of M — but for most of this paper, these will not play an important role. 
Since we are interested in topology change, we shall use path integral techniques to 



quantize the action (|2.1| ). But because we are dealing with manifolds with boundary, we 
must first determine the appropriate boundary conditions. A simple heuristic argument is 
as follows. Let us start by picking boundary values for either e or uj on dM. This amounts 
to choosing Dirichlet boundary conditions for (say) uj. Now, the kinetic term in the action 
(|2TT| ) can be written as < e, *duj >, where * is the Hodge star operator and <, > is the 
corresponding inner product. But if uj obeys Dirichlet boundary conditions, *duj obeys 
Neumann boundary conditions, so we ought to require the same of e in the inner product. 
We thus expect opposite boundary conditions for our two fields. 

To make this argument more rigorous, let us first examine the Hilbert space of (2+1)- 
dimensional quantum gravity. Canonical quantization of this theory on a manifold RxE 
has been discussed by a number of authors; see |7j for a summary. A key feature is that the 
classical reduced phase space is a cotangent bundle, whose base space is the space of flat 
connections wona slice S modulo SO (2,1) gauge transformations and large diffeomorphisms. 
In the simplest approach to quantization, the "connection representation" ||, states are 
therefore gauge-invariant functionals \l/[u;j a ] of the spatial part of the connection, subject to 
the constraint that uj be flat on E. 

The corresponding boundary conditions for the path integral therefore require us to to 
fix a flat connection uji a on dM. More precisely, let /: dM^-M be the inclusion map. We 
can then freely specify the induced connection one-form I*uj on dM, as long as the induced 
curvature I*R vanishes. SO(2,l) gauge invariance of the resulting amplitude is formally 
guaranteed by the functional integral over the normal component of uj: at dM, uj± a is a 
Lagrange multiplier for the constraint 

N a = l e V (d iej a - d jei a + e abc (uj lb e JC - u ic e jb )) (2.6) 

that generates SO (2,1) transformations of I*u ||||. Observe that we must integrate over 
uj± at the boundary to enforce this constraint — that is, we must not fix uj± as part of the 
boundary data. This is in accord with the canonical theory, in which wave functionals depend 
only on the tangential components of uj. 

Note that the flat connection I*uj is completely determined by its holonomies, that is, by 
a group homomorphism 

H : Tii(dM) — ► SO(2,l). (2.7) 
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As we shall see later, the transition amplitude can be described rather explicitly in terms of 
these holonomies. 

The specification of I*uo is not quite sufficient to give us a well-defined path integral. To 
obtain an additional boundary condition, it is useful to decompose u into a background field 
uj that satisfies the classical field equations and a fluctuation Q: 

u) = ou + n, duo a + ^ abc ^b Aw c = 0. (2.8) 

Note that if there is no classical field with our chosen boundary values, then the transition 
amplitude is zero: the integral over e gives a delta functional 5[R a ] that vanishes everywhere. 
The requirement of existence of a classical solution leads to Witten's selection rules for 
topology change |4[]. Assuming now that Co exists, the boundary condition I*Q = can 
be recognized as a part of the standard Dirichlet, or relative, boundary conditions for a 



one-form 10,11 



Ptt = = *D*Sl. (2.9) 

Here, * is the Hodge star operator with respect to an auxiliary Riemannian metric h that 
we introduce in order to define a direction normal to the boundary, while D is the covariant 
exterior derivative coupled to the background connection u, 

Df3 a = d(3 a + e abc u b A (3 C . 

These boundary conditions make the operators D and *D* hermitian conjugates, and guar- 
antee that the Laplacian A = D*D* + *D*D is hermitian. We shall see below that the new 
condition on the derivatives of Q, which can be written in component form as D 11 VI^\qm = 0, 
is actually a gauge condition. Note that fl2.9p depends on the nonphysical metric h; we must 
check that the final transition amplitudes are independent of h. 

We next turn to the boundary conditions for the triad e. I*e and I*u are conjugate 
variables, so we cannot expect to specify them simultaneously. On the other hand, we must 
not integrate over the normal component of e at the boundary. Indeed, e± a acts as a Lagrange 
multiplier for the constraint 

A> = ^e ij (d t u 3 a - djUi 11 + e abc u ib u JC ) , (2.10) 

and would lead to a delta functional 5[N a ] = 5[I*R a ] at the boundary. But we have already 
required that I*u be flat, so such a delta functional would diverge. We avoid this redundancy 
by fixing e± at DM. We shall see below that transition amplitudes do not depend on the 
specific value of e±, so this does not contradict the canonical picture. 

As with uj, we can obtain additional boundary conditions by decomposing e into a classical 
background field and a fluctuation 

e = e + E, de a + e abc to b A e c = (2.11) 

where E± vanishes, i.e., I*(*E) = 0. This restriction on E is a part of the standard Neumann, 
or absolute, boundary conditions for a one-form, 

I*(*E) =0 = I*(*DE). (2.12) 

Once again, these conditions make the Laplacian A hermitian. 



4 



3. Path Integrals and Ray-Singer TorsionT] 



We are interested in path integrals of the form 

Z m [Flu} = J [duj}[de] exp {*/ grav [M]} , (3.1) 

where M is a manifold whose boundary 

<9M = £iH£ 2 (3.2) 

is the disjoint union of an "initial" surface Si and a "final" surface S 2 . (Si and S 2 need not 
be connected surfaces.) Our first step is to choose gauge conditions to fix the transformations 



( |2. 1| )-( |2~5| ). To do so, we employ the auxiliary Riemannian metric h introduced in the last 
section, and impose the Lorentz gauge conditions 

*D*E a = *D*Q a = 0. (3.3) 

For later convenience, we use the covariant derivative D coupled to the full connection u 
rather than u in our gauge-fixing condition. D and D agree at the boundary, however, so 
the gauge condition on Q reduces to the second equation of fl2.9|) on dM, as promised. 

To impose (|3.3| ) in the path integral, we introduce a pair of three-form Lagrange multi- 
pliers u a and v a , and add a term 

W = - / («« A *D*E a + v a A *D*Q a ) (3.4) 

JM 

to the action. It is not hard to see that for the path integral to be well-defined, u should 
obey relative boundary conditions (I*(*D*u) = 0), while v should obey absolute boundary 
conditions (*t> = on dM). The latter restriction again has a rather straightforward inter- 
pretation: since we are already imposing the gauge condition (|2.9|) on Q at the boundary, 
we do not need the added delta functional <5[*.D*f2] that would come from integrating over 
v at dM. 

As usual, the process of gauge-fixing leads to a Faddeev-Popov determinant, which can 
be incorporated by adding a ghost term 

W = - J (/A *D*Df + gA *D*Dg) , (3.5) 

where /, /, g, and g are anticommuting ghost fields. We must again be careful about 
boundary conditions: corresponding to restrictions and ( |2.12|) on Q and E, we choose 
/ and / to satisfy relative boundary conditions and g and g to satisfy absolute boundary 
conditions. The full gauge-fixed action is then 

-^grav -^gauge -^ghost 

= J [E a a (Dn a + x -t ahc n b Att c + *D*u a 

+ ^e abc e a AQ b AQ C -v a A*D*Q a - f A*D*Df -g A*D*Dg]. (3.6) 



*Much of this section is a summary and elaboration of Witten's work in 
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E and v occur linearly in Q3.6|) , so following Witten [|J, we shall first integrate over these 
fields to obtain delta functionals. There is one subtlety here: certain modes of E do not 
contribute to the action, and must be treated separately in the integration measure. The 
issue is most easily understood in the case of a linear action, for instance 

M 

If we expand a and (3 in terms of orthonormal modes <p n of the Laplacian, a = Y^Un^m 
(3 = J2b n (pn, it is easy to see that / = X/A n a n & n , where A n are the eigenvalues of the 
Laplacian and the sum automatically excludes the coefficients of the zero- modes of A. The 
path integral measure [da] thus splits into an integral over modes for which A n 7^ — giving 
a delta function n'^n 1( ^(^n) — and an integral /[<iao] over zero-modes. 

Unfortunately, the action ( |3.6| ) is not linear in Q, and such a mode decomposition is much 
more difficult.]] However, we can argue as follows. The integral over the "nonzero-modes" of 
E will give a delta functional of DQ a + \e abc VL b A Q c + *D*u a . The zeros of this expression 
form a surface (Q(s),u(s)) in the space of fields, and if we expand the action around these 
zeros, only those fields infinitesimally close to this surface should contribute to the path 
integral. Writing Q = Q + 5Q, we easily find that the relevant zero-modes of E are those E 
for which 

D Q+Cl E = 0. (3.7) 

(Note that the E depend on Q, so the order of integration below cannot be changed.) 
Performing the integration over E and v, we obtain 

[dQ}[du][dE}[dv}e u = J [dn}[du][dE]5[DQ a + ^e abc Q b AQ C + *D*u a ]5[*D*Q a }. (3.8) 

The argument of the first delta functional vanishes only when D*D*u a = 0; assuming that 
the connection u is irreducible, this implies that u a = 0. The delta functional then imposes 
the condition 

Dn a + ^e abc n b Att c = o, (3.9) 

which can be recognized as the requirement that u = Co + Q be a flat connection. This, in 
turn, allows us to eliminate the term 

/ -tabc e a A Q b A Q c = - I e a A DQ a = f De a AQ a = 
Jm 2 Jm Jm 

in®. 

We can now use the delta functionals to perform the remaining integration over Q. It is 
straightforward to show that 



[dn]s 



DQ a + ^e abc Q b A Q c + *D*u a 



5[*D*Q a ] = [du)]|det'i7 el |-\ (3.10) 



^See 113] for a related treatment of the Abelian B-F system, in which this nonlinearity is not an issue. 
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where uj = Co + Q ranges over flat connections with our specified boundary values and 
L rcl = *Dq + D&* maps a one-form plus a three-form (a, (3) obeying relative boundary 
conditions to a one-form plus a three-form (*D^a+D^*P, D^*a) obeying absolute boundary 
conditions. Performing the ghost integrals, we finally obtain 

r _ det'A^det'A^ 8 

z m [im = J mm ® (o) , (3.H) 

where A(^) is the Laplacian *D^*Dcj + D^*D^j* acting on fc-forms. 

(Strictly speaking, the determinants det'A are not well-defined for a noncompact gauge 



group like SO (2,1). But this is a minor problem, whose solution has been discussed in [O 



and |L4| , and it does not affect our final expression for amplitudes in terms of Reidemeister 



torsion |T3| .) 

Now, by expanding one- forms and three-forms in modes of LLL_, one may easily show 
that 

|(det'L_)(det , L t _)| = det'Agdet'Ag = det'A^det'A^ 8 (3.12) 

Moreover, det'A^ = det'A^^, since the Hodge star operator maps any eigenfunction a of 

A/|\ to an eigenfunction *a of Agf.^ with the same eigenvalue. Some simple manipulation 
then shows that 



Z M [I*u] = J[dCb][dE]T[Cd] 



(det'Ag^det'Ag)^ = (det-Agp^^det^A^) 1 / 2 ; ^ } ^ 

(det'Ag) " (det'Ag- 



The combination of determinants T[u] may be recognized as the Ray-Singer torsion [Kj; the 



equality of relative and absolute torsions on odd- dimensional manifolds is shown in ||15|| . A 
similar transition amplitude for Abelian B-F theories has bee discussed by Wu [IB]. 



We can now return to the question of whether our amplitude Zm depends on the choice 
of auxiliary metric h. When there are no zero-modes, T[u] is known to be independent of 
h |]ID| , |r7j . When zero-modes are present, T[u] depends on h, but so does the volume element 
[du] [dE] . If, as we have assumed, there are no ghost zero-modes — that is, if H°(M;V^) = 
H°(M,dM; V a ) = 0— then the combination [duj}[dE}T[ud} is again independent of h [0,[19|. 
If ghost zero-modes are present, they must be included in the integral (|3.13| ); when they are, 
the amplitude is once again independent of h. 



We conclude this section with a discussion of the range of integration in ( |3.13j ). As noted 
above, u ranges over the space of gauge-fixed flat SO(2,l) connections — i.e., the moduli 
space of flat connections modulo gauge transformations — with specified boundary values. In 
general, this space is rather complicated!] We can at least determine its dimension, however, 

■^For the remainder of this paper, we shall treat this space as if it were a manifold. When M has the 
topology RxE, this is essentially correct pOfl ; for more complicated topologies, we do not know whether this 
assumption is justified. 
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by linearizing ( fT3| ) and (|3.9|) : solutions <5f2 of 

D Q 5tt a = = * 5Vl a (3.14) 



satisfying the boundary conditions ( |2.9|) are cotangent vectors to the moduli space of flat 
connections at uj. 

These conditions have a natural cohomological interpretation. Since the connection uj 
is flat, D\ = 0, so we can construct a de Rham cohomology H*(M; V&) on the complex of 
forms on M with values in the flat bundle determined by uj. Then ( |3.14| ) is equivalent to 
the condition 

5tt E iif 1 (M, dM; (3.15) 

where the boundary conditions determine the use of relative cohomology. The range of E 
has a similar interpretation. From ( |3.3j ) and (|3.7| ) and the boundary conditions (|2.12| ), it 
follows that 

E E if X (M; Vq). (3.16) 

Because of the change in boundary conditions, E and 5Q do not lie in the same cohomology 
groups; indeed, by Poincare-Lefschetz duality, 

H\M- V a ) « H 2 (M, dM; V a ). (3.17) 



Unlike (|3.15 ), equation (|3.16| ) determines the entire space of fields E, not merely its 



tangent space. Moreover, the integrand T[uj] in ( p.!3| ) is independent of E. This means that 
if the cohomology group ( |3.17| ) is nontrivial, the amplitude Zm[I*w] diverges. This is the 
infrared divergence cited by Witten as an indication of classical behavior [|J. In contrast to 
the cases discussed by Witten, however, our boundary conditions allow this divergence to 
appear even when the moduli space of flat connections consists of isolated points. This 
will be the case in specific examples we discuss below. 



4. From Ray- Singer to Reidemeister Torsion 



In principle, the integral (|3.13| ) determines the transition amplitude for an arbitrary 



topology change in 2+1 dimensions. In practice, however, the determinants in T[uj] are 
usually impossible to evaluate. We must therefore take one further step, and relate the 
Ray-Singer torsion to the combinatorial, or Reidemeister, torsion. 

We begin with a brief description of the Reidemeister torsion. (For more details, see JTO 



plf , or 1 2~2|]). It is instructive to start with a concrete description of the flat bundle Let 



7Ti be the fundamental group of M, and let M denote the universal covering space of M, so 

M»M/jti. (4.1) 

As in ( |2.7|) , the flat connection uj is determined up to gauge transformations by its holonomy 
group 

H : 7Ti -> SO(2,l). (4.2) 
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We can now define our flat bundle as 



V a = (Mx 80(2,1))/^, (4.3) 



where 71*1 acts on M as in (4.1) and on the Lie algebra so(2, 1) by the adjoint action of the 
holonomy group. 

We now repeat this construction in a slightly different form. Let us treat M as a CW 
complex, with fc-cells {e? k ^}. M can be viewed as a fundamental domain embedded in M, 
which then has a corresponding cell decomposition in terms of the translates {ge? k ^, gEiti}. 

The chain groups Ck(M) of the universal covering space M thus become modules over the 
real group algebra R[7Ti], with the constituting a preferred basis.f] Note that if is 

a cell in M, its boundary will not, in general, lie in M, but will rather be a sum of translates 
of (k — l)-cells by elements of 7i"i. Relative to our preferred basis, the boundary operator 
can thus be viewed as a matrix with elements in H[7Ti]. 
We can now construct the twisted chain complex 

C(M; V a ) = C(M) ® R[ltl] so(2, 1), (4.4) 

where, as in Q4.3| ), ~K\ acts on C(M) by translation and on so(2, 1) by the adjoint action. 
A preferred basis for C(M; Vq) consists of the elements {e" fe ) <8> J a }, where the J a are an 
orthonormal set of generators of so(2, 1). We shall abbreviate these basis elements by c^) 
below. The boundary operator for this twisted complex can be viewed as a real matrix: if 

de (k) = 0effc-i) + ■ • •> then 

5(4) ® J a ) = ef fc _ 1} ® ac%)J a + . . . = ef fc _ 1} ® 5[^] a fe J fe + . . . , (4.5) 

where 5'[(7] a fe are the matrix elements of g in the adjoint representation. 

The Reidemeister torsion of M is now defined as follows. We have chosen a preferred basis 
C(k) for each of the chain groups Ck(M; V&). We also have a preferred basis for each homology 
group i?fe(M;V&), determined from the harmonic forms of the previous section by the de 
Rham isomorphism. Let us select a set h/w £ to represent these basis elements. We next 
choose an arbitrary basis b(fc) for each B k , and a set b^-i) in Ck such that db^-i) = t>(fe— i)- 

It is easy to see that the set (b( fc ), h(k), k(&-i)) forms a new basis — call it £( fc ) — for Ck- 
Let T(fc) denote the matrix representing the change of basis from c^) to c/m, that is, C(^) = 
T(k) c {k)- The Reidemeister torsion is then defined as 

det Tcoidet T/ 2 \ 

t(M;Vz) = (0) (2) . 4.6 
detT(!)detT( 3 ) 

The relative Reidemeister torsion r(M, 9M; V^) is obtained by the same construction with 
the relative chain complex. Since cellular decompositions of manifolds are often rather 
simple, this invariant can sometimes be calculated quite explicitly. 

*As usual, we denote chain groups by Ck, the kernel of d in Ck by Zk, and the image of Ck+i in Zk by Bk] 
the homology groups are Hk — Zk/Bk- 
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For a manifold without boundary, the Reidemeister torsion is, remarkably, equal to the 
Ray-Singer torsion T[u] of equation ( ft.l^) P2| , P3[ . When a boundary is present, this equality 



no longer holds, but only a small correction is needed fi5| : 

T[u>] = 2 3 * (aM)/4 r(M, 8M- V a ), (4.7) 

where x(9M) is the Euler characteristic of the boundary (The factor of three in the exponent 
is the dimension of so(2,l).) We can therefore determine the transition amplitude ( |3.13|) by 
evaluating the Reidemeister torsion for M. 



5. Topology- Changing Manifolds 

To compute a topology-changing amplitude, it remains for us to specify the manifold 
M that mediates the transition. This choice is not trivial — as Amano and Higuchi have 
shown [||], the requirement that the initial and final boundaries be spacelike strongly restricts 
the allowed topologies. We begin with M 2 , the simplest of the interpolating manifolds of 
reference ||, which describes a transition from a genus three surface to two genus two 
surfaces. This manifold, shown in figure 1, can be described as follows: 

Let V be a genus four handle-body, and let V denote a reflected copy of V. Remove 
from the interior of V two genus two handle-bodies, as indicated by the shaded areas in 
figure 1, to obtain a manifold W\ with boundary dW\ ~ dVU S 2 U S' 2 , where S 2 and S' 2 are 
genus two surfaces. Next, remove from the interior of V a single genus three handle-body, 
as shown, to obtain a manifold W% with boundary dW 2 ~ dV U S 3 , where S 3 is a genus 
three surface. Now identify W\ and W 2 along their common boundary dV ~ dV to obtain 
our desired manifold M 2 , which clearly has a boundary dM 2 ~ S 2 U S 2 U S3. More general 
amplitudes can be obtained by combining manifolds M n that mediate between a genus n + 1 
surface and n genus two surfaces; we refer the reader to reference || for details. 

Figure 1 also indicates a set of generators pi, . . . , p 8 for 7r 1 (M 2 ), which obey the relations 



[P2,Pl 1 ][P 3 1 >P4 1 ] = [P 6) P5 1 ][P7 1 ) P8 1 ] = !> 

P3=P5, (5.1) 

where [a, t] = oto~ x t~ x is the commutator. Apart from a slightly unusual normalization 
that will be useful later in describing the cell decomposition, the relation obeyed by pi, . . . , P4 
can be recognized as that of a Fuchsian group that uniformizes a genus two surface |24j . The 
generators p 5 , . . . , p 8 determine a similar group. 

It will be helpful to compare the amplitude coming from M 2 to a related topology- 
preserving amplitude. We obtain the latter from a new manifold P 2 constructed by attaching 
two copies of W2 along their common boundary dV . This manifold mediates transitions 
from S 3 to another genus three surface S' 3 . Figure 2 shows P 2 , along with a set of generators 
01, . . . , <r 8 of ni(P 2 ), which satisfy the relations 

^J^oT^kiTVl V4 1 = ^eVs Vf^t^Veke = °5- ( 5 - 2 ) 
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To compute the Reidemeister torsion, we need a cell decomposition for the universal 
covering space M 2 . The first step in obtaining such a decomposition is to dissect M 2 into 
a simply connected region Um that can serve as a fundamental domain upon which iti(M 2 ) 
acts by translation. We do so by cutting M 2 open along a set of surfaces transverse to the 
generators of the fundamental group. 

Since this process is somewhat difficult to visualize, let us begin with a simple example. 
Let Ej i be a genus one surface with a single hole, and consider a thickening of E^i into a 
three-manifold N, one of the elementary building blocks of M 2 . Figure 3 shows a dissection 
of N into a simply connected fundamental domain. It is clear that N can be recovered by 
identifying the boundaries of this figure; this identification represents the action of irx(N). 

Figure 4 shows the corresponding dissection of M 2 . The cuts on the top half of the figure 
are evidently analogous to those of figure 3; they are then extended through the boundary 
dV into the bottom half of the figure, where they are continued until they reach the inner 
boundary E 3 . 

The corresponding construction for P 2 is illustrated in figure 5. We start with two 
copies of the bottom half of figure 4. However, a new noncontractible loop now appears, 
representing the element 05 of tti(P 2 ), that was not present in figure 4. To obtain a simply 
connected fundamental domain Up for P 2 , we therefore need one more cut transverse to this 
loop. This final cut, required by the larger fundamental group of P 2 , will play a critical role 
in the suppression of topology change — it will lead to the existence of extra zero-modes, and 
thus more highly divergent infrared behavior, in the topology-preserving amplitude. 

To obtain cell decompositions of M 2 and P 2 , it now suffices to find cell decompositions of 
U m and Up compatible with the action of their fundamental groups — that is, decompositions 
for which the process of gluing (say) Um back together to obtain M 2 identifies like cells. 
The final results are shown in figures 4 and 5. We have labeled a basis e%\ of fc-cells in 
each figure; the remaining cells are translates of these basis elements by elements of the 
appropriate fundamental group. 

6. An Explicit Computation 

We are finally ready to calculate the amplitude Zm 2 for tunneling from the genus three 
surface E 3 to the genus two surfaces E 2 and E' 2 . We shall proceed in two steps, first obtaining 
some general information about the uj and E integrals in ( |3.13| ) and then computing the 
Reidemeister torsion for a particular choice of flat connection. 

We begin by showing that the topology-preserving path integral Zp 2 has at least six 
more zero- modes than Zm 2 - This is in itself enough to indicate a suppression of topology 
change — each E mode leads to an infrared divergence, so this mismatch shows thatZp 2 is 
more divergent than Zm 2 - We caution the reader, however, that the choice of P 2 to "normal- 
ize" M 2 is somewhat arbitrary; until we understand more about the overall normalization 
of amplitudes, we can draw no firm conclusions about absolute probabilities for topology 
change. 
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Observe first that the twisted Euler characteristics p5 | 

J2(-iy^mH t (M 2 , dM 2 ; V a ) = *(M 2 , dM 2 ; V a ) = ^ X (dM 2 ; V Q ) = 12 

£(-l)VM,(P 2 , dP 2 ; V a ) = X (P2, dP 2 ; V a ) = \x(dP 2 ; V ) = 12 (6.1) 

serve as constraints on the dimensions of the twisted homology groups, placing a lower limit 
of twelve on the total number of zero- modes. In the case of M 2 , this limit is realized. Indeed, 
the twisted cell complex associated with M 2 is 

C (M 2 ,dM 2 ;V„) & diM^dM^Vc) & C 2 (M 2 , dM 2 ; V a ) & C 3 (M 2 ,dM 2 ;V„). 
£- R©...©R R©...©R £■ R © R © R 

V v ' " v ' 

6 times 21 times 

(6.2) 

If we can show that d 2 is an epimorphism and that d 3 is a monomorphism, it will fol- 
low that H (M 2 ,dM 2 ;V a ) = H 1 (M 2 ,dM 2 ;V s> ) = H 3 (M 2 ,dM 2 ;V Q ) = 0, and hence that 
dimH 2 (M 2 ,dM 2 ;Vz) = 12. 

Given a set of generators { J°, J 1 , J 2 } for SO(2,l), the adjoint action described in section 
4 determines a matrix representation of 7Ti(M 2 ), 

S: 7n(M 2 ) -> GL(3,R). 

To show that d 2 is an epimorphism, consider the following six boundary operations, which 
can be read off figure 4: 

9 2 ( e ; 2) ©j a ) = e ; 1) ©(i-% 3 ])j a 

d 2 {e\ 2) © J a ) = ef 1} © % 2 ](1 - S[ Pl })S[p 3 VrV 3 K a - (6-3) 

The S[pi\ are nondegenerate matrices which each stabilize a one- dimensional subspace of R 3 , 
so 1 — S[p3] is a rank two matrix that is zero on the subspace stabilized by S[ps\. Similarly, 
^[/^(l — S[px])S[p3 x p~l 1 p 3 ] is a rank two matrix that is zero on the subspace stabilized 
by S' = S[p3 V4P3P1P3 1 Pi Va]- But in order for the boundary dM 2 to be spacelike, the 
commutator [p 3 , p% V4P3P1P3 1 p1 1 p 3 ], which represents a loop on dM 2 , must be nontrivial. 
Hence S[p 3 ] and S' stabilize different one-dimensional subspaces, and the images of 1 — S[p 3 ] 
and 1 — 5" span all of R 3 . This, in turn, means that the image under d 2 of the space spanned 
by {eLj © J a , e 3 2 ) © J a } contains {eh\ © so(2, 1)}. A similar argument shows that the image 
of the space spanned by {e%-> © J a , ejL © J a } contains {e 2 ^ ©so(2, 1)}, thus establishing that 
d 2 is an epimorphism. 

An analogous argument, starting with the boundary operations 

ds(e { 3) © J a ) = e\ 2) © (1 - S\pl ']) J a - e\ 2) © (1 - S\f%\)J a + ... (6.4) 

and the nontriviality of the commutator [pj , p 3 ], shows that 83 is a monomorphism. We 
have thus established that H 2 (M 2 , dM 2 ; Vq) is twelve-dimensional. This conclusion has been 



12 



checked by making explicit choices for the holonomies {^[pj]} and identifying generators of 
the twisted homology groups. 

We can repeat the same analysis for P 2 , but let us instead take a shortcut. It is easy to 
see that H (P 2 , dP 2 ; Vq) = 0, so the constraint (|Q| ) implies that 

dimH 2 (P 2 , dP 2 ; Va) > 12 + d%mE x {P 2 , dP 2 ; V a ). 

So if we can show the existence of three generators of Hi(P 2 , dP 2 ; Vq), we will have proven 
that Zp 2 has at least six more zero-mode integrations — three Q modes and three E modes — 
than Z M2 . 

But this is apparent from figure 5: no combination of (g> J a and e%\ (g> J a lies in the 
image of d 2 , and 

di (e 3 (1) ®J a + 4) ® J a ) = 0. (6.5) 

Hence + ef^) (g) J a } represent three generators of Hi(P 2 , dP 2 ; Vq), as claimed. 

We conclude this section by describing an explicit computation of the Reidemeister tor- 
sion for M 2 . This requires a choice of flat connection uj, which may be determined by its 
holonomies around the curves pi that generate Ti\{M 2 ). Note, however, that each of the pi 
can be deformed to a curve on one of the boundary components of M 2 . The connection uj 
is therefore fixed by its holonomies on S 2 U E 2 U S 3 , and hence by the boundary data Puo. 
This means that the relevant moduli space of flat connections is a single point, and that the 
integral over uj disappears from ( |3.13|) , as expected from the vanishing of H 1 (M 2 , dM 2 ; 



Following reference |J, we choose the S[pi] as follows. S[pi] through S\p^] are the gener- 
ators of an arbitrary Fuchsian group uniformizing a genus two surface. Such groups form a 
six-parameter family, which can be written down from Fenchel-Nielsen coordinates by using 
the results of [p6| . A convenient two-parameter family S[pi](k, r) is given in the appendix. 
We then take 

% 5 ] = % 3 ]> 5 [P 6 ] = % 4 ]> 

s [pj] = 5[p 4 P 2 PiP 2 1 Pi 1 }, S[p 8 ] = S[p A p 2 p± l }. (6.6) 

This expression differs slightly from that of reference 0, but only because of our different 
choice of generators of 7r 1 (M 2 ). Amano and Higuchi show that with these choices, the 
boundaries of M 2 are spacelike and nonsingular. 

As described in section 4, calculating the torsion requires computing the determinants of 
the matrices T^) that give the change of basis from cn.) to (bo,), h(fc), b(fc_i)). The total mea- 
sure appearing in the amplitude ( |3.13| ) is independent of the choice of basis (b(k), h(fe)j b(k-i)), 



but the determinants ( [4.6|) by themselves are not — we must make an explicit choice of the 
homology basis ha.) for a numerical value of the torsion r to have meaning. For now, we 
choose the simple but rather arbitrary basis h^) described in the appendix. A tedious but 
routine calculation then gives a Reidemeister torsion as a complicated rational function of 
k and r. A plot is shown in figure 6. The torsion falls off very rapidly for large values of k 
and r; for example, 

r ( r = io, fc = 1000) «7x 10" 52 . 
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To interpret these results, we must understand the remaining E integrals in ( [3.131 ). The 
measure dE is determined by the homology basis h( 2 ). We first select an orthonormal basis 
E a such that 




where is the union of one-cells dual to and the integral includes a trace over the 
generators of so(2,l) (see |T(| for details). If we then decompose an arbitrary zero-mode as 
E = c a E a , the integral ( |3.13| ) is / dc a , and the dependence of c a on the choice of homology 
basis cancels the dependence of the torsion r. 

If we wish to cut off the infrared divergent E integrals, however, this basis dependence 
reappears — the range of integration will, in general, depend on the choice of h( 2 ). This 
dependence can be translated into a dependence on the auxiliary metric h introduced in 
section 3 to fix the gauge. The difficulty appears to be a typical regularization problem — we 
do not know how to regulate our integrals in a way that respects the invariance of the original 
theory. 

There is some hope for a physical resolution of this problem. Witten has suggested that 
the infrared divergences of transition amplitudes in (2+l)-dimensional gravity reflect the 
appearance of infinite- volume "classical" spacetimes. If we could give a concrete geometrical 
meaning to the limits of integration, cutting off at some (observable) scale, we might be able 
to define a basis-independent regularization. It is also possible that the addition of matter 
to our vacuum theory might regulate the divergences, again by limiting maximum lengths. 
We leave these questions for future investigation. 

7. Conclusion 

We have now seen that path integrals representing spatial topology change in (2+1)- 
dimensional general relativity need not vanish. Starting with any initial data in our two- 
parameter family — or, by a long but straightforward generalization, any other admissible 
initial geometry — we can explicitly compute the torsion T[I*u], and thus the amplitude 
( |3.13| ). Indeed, we have seen that such topology-changing amplitudes may diverge, thanks 
to the existence of zero- modes E a of the triad e a . These divergences presumably reflect the 
appearance of "classical" spacetimes, in which distances measured with the metric g^ v = 
Gfj, a e ua become arbitrarily large 0. 

Nevertheless, our results may be interpreted as providing evidence that topology change 
is suppressed. We have seen that while the topology-changing amplitude mediated by M 2 is 
infrared divergent, the closely related topology-preserving amplitude mediated by P 2 is even 
more divergent. Evidently, no firm conclusions about topology change can be drawn without 
a much better understanding of the overall normalization of amplitudes in (2+l)-dimensional 
gravity. This is a difficult problem: we must not only consider an infinite number of possible 
interpolating manifolds, but must also find a sensible way to regulate infrared divergences 
without breaking the symmetries of the original theory. Clearly, much work remains to be 
done. 
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One important step would be to find an easier and more intuitive method for computing 
the degree of divergence for an arbitrary interpolating manifold without requiring a full cell 
decomposition. We do not have a complete answer to this problem, but we offer the fol- 
lowing observations. As we have seen, Hi(M,dM; V u ) counts the dimension of the moduli 
space of flat connections on M with specified boundary data. In general, a flat connection 
is determined up to gauge transformations by its holonomies, that is, by an assignment of 
an element of SO (2,1) to each independent generator of 71*1 (M). For our manifold M 2 , it 
is evident from figure 1 that every generator of the fundamental group can be deformed 
to the boundary; thus, the connection is completely determined by boundary data, and 
dimHi(M2,dM2;V UJ ) = 0. For P2, on the other hand, one independent generator — 04, for 
example — cannot be deformed to the boundary, and accounts for our three generators of 
ifi(M 2 , dM 2 ; VJ). We do not know whether this argument can be made rigorous when gen- 
eralized to an arbitrary three-manifold, but such an extension may be possible. Similarly, 
it may be possible to obtain information about dimH2(M,dM;V tJ ) = dimH\{M;Vu) by 
looking at connections on the double of M. Note that paths connecting boundary com- 
ponents of M become closed loops on the double, and may contribute to the dimensions 
dimH*(M,dM;Vu)- Finally, the Euler characteristic constraint ( |6.1|) holds for any three- 
manifold, and places a useful constraint on the number of divergent integrals. 
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Appendix. Computational Details 

In this appendix, we briefly describe some of the intermediate steps in the computation 
of the torsion plotted in figure 6. We begin with a choice of holonomies, which we give in 
the two-dimensional representation of SL(2,1R): 



These describe a genus two surface consisting of two identical one-holed tori (£1,1 of figure 
3) joined with a relative twist; r parametrizes the twist, while k parametrizes the length of 
a closed geodesic on each copy of £1,1. 




(A.l) 
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We next describe the choice of homology basis used in the computation of the torsion 
r. Using the cell structure shown in Figure 4, we begin by selecting basis elements b(i) and 
b( 2 ) as defined in section 4: 

4) = e (2) ® J°i "(1) = e (2) ® J°i 4) = e (2) ® 

%j = ef 2) ® J°, % )= eJ 2) ®J°, bf 1} = e\ 2) <g> J°, (A.2) 

~ b f 2) = e (3) ® J a , a = 1,2,3. (A.3) 

Our homology basis is then 



ty 2 ) 


= 4) 


® J 2 + 


«(2) 


= 4) 5 


5 ^+5:^1). 


1(2) 


= 4) 




1(2) 


= 4) s 




ft 5 


- e (2) 


® j 1 + E^w» 


ft 6 

1(2) 


= 4) ^ 




1(2) 


= 4) 


® J 1 + 


1(2) 


= 4) $ 


$J 2 +E^(D> 
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= 4) 
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1(2) 


= 4) 


§j 2 + E^(d 


ft 11 - 

1(2) " 


= 4) ( 




h 12 

1(2) 


= 4) 5 


5J 2 + E^(i) 



(A.4) 

where the coefficients K l a , whose exact values are not needed for the computation of the 
torsion r, are uniquely determined by the requirement that dh? 2 \ = 0. 

Of course, this choice of homology basis depends on the boundary operator d, and there- 
fore on the holonomies pj. The choice ( |A.4| ) appears to be valid for generic values of the 
holonomies, but there are points at which linear combinations of our h%) lie in the image 
of 83. At these points, (|A.4|) is no longer a valid basis, and our computed value of the 



Reidemeister torsion r is not correct — in fact, it appears to go to zero. This behavior is an 
artifact of our basis choice, and does not affect the integral (|3 . 1 3|) . 



We do not include the final functional form of r(r, k) — it would require three pages — but 
it is available from the authors. 
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Figure 1. The manifold M2 is formed by identifying W\ and W2 along their common genus four 
boundary. A basis of loops for tt^M^) is shown. 
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Figure 3. A dissection of the thickened one-holed torus into a simply connected fundamental 
domain. 
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Figure 4. A dissection of M^. A basis of cells efy is shown. 
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Figure 5. A dissection of P2 and a basis of cells. 
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Figure 6. A two-parameter family of Reidemeister torsions r(r, k) for different choices of the flat 
connection on dM2', the range shown is 1.1 < r < 12, .3 < k < 1.2. 




